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ABSTRACT 

A new class of X-ray binaries has been recently discovered by the high energy observatory, IN- 
TEGRAL. It is composed of intrinsically obscured supergiant high mass X-ray binaries, unveiled by 
means of multi-wavelength X-ray, optical, near- and mid-infrared observations, in particular photo- 
metric and spectroscopic observations using ESO facilities. However the fundamental questions about 
these intriguing sources, namely their formation, evolution, and the nature of their environment, are 
still unsolved. Among them, ICR J16318— 4848 - a compact object orbiting around a supergiant B[e] 
star - seems to be one of the most extraordinary celestial sources of our Galaxy. We present here new 
ESO/VLT VISIR mid-infrared (MIR) spectroscopic observations of this sourc^E First, line diagnos- 
tics allow us to confirm the presence of absorbing material (dust and cold gas) enshrouding the whole 
binary system, and to characterise the nature of this material. Second, by fitting broadband near 
to mid-infrared Spectral Energy Distribution - including ESO NTT/Sofl, VLT/VISIR and Spitzer 
data - with a phenomenological model for sgB[e] stars, we show that the star is surrounded by an 
irradiated rim heated to a temperature of ~ 3800 — 5500 K, along with a viscous disk component at 
an inner temperature of ~ 750 K. VISIR data allow us to exclude the spherical geometry for the dust 
component. This detailed study will allow us in the future to get better constraints on the formation 
and evolution of such rare and short-living high mass X-ray binary systems in our Galaxy. 

Subject headings: circumstellar matter — infrared: stars — stars: emission-line, B[e] — X-rays: 
binaries — X-rays: individual (IGR J16318-4848) 



1. INTRODUCTION 

The high energy INTEGRAL (INTErnational Gamma- 
Ray Astrophysics Laboratory) observatory, observing in 
the 20keV— 8MeV range, has performed a detailed sur- 
vey of the Galactic plane, and the ISGRI detector on the 
IBIS imager has discovered many new sources, most of 
them reported in iBird et al.l ([2010ft . Many of them are 
concentrated in a dire ction tangent to the Norma arm 
of th e Galaxy (see e.g. iTomsick et al.l [20041 : IChatv et al.l 
120081) . rich in star forming regions. The most impor- 
tant result of the INTEGRAL observatory to date is 
the discovery of many new high energy sources exhibit- 
ing common characteristics which, previously, had rarely 
been seen. Most of them are High Mass X-ray Binaries 
(HMXBs) hosting a neutron star orbiting a supergiant 
O/B companion star. Some of the new sources are very 
obscure d, exhibiting huge intrinsic and local extinction 
(see e.g. lChatvll2011ft . Certainly one of the most extraor- 
dinary examples among the obscured high energy sources 
is the extremely absorb ed source IGR J16318— 4848 (see 
IFilliatre fc CIiatvl[200l . 

IGR J16318— 4848 was the first source to be discovered 
by ISGRI on 29 January 2003 at the Galactic coordinates 
(/,&) ~ (336°, 0.5°) with an uncertainty radius of locali- 
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sation of 2' jC ourvoisier et all 120031 ). ToO observations 
were then triggered with XMM-Newton, allowing a more 
accura te localisation (RA =16:31:48.6; Dec=-48:49:00) 
at 4" ([Walter et al.l 12003ft . XMM-Newton observations 
showed that the source was ex hibiting a strong column 
density of A H ~ 2 x 10 24 cm" 2 ([Matt fc Guainazzill2003l : 
I Walter et al"]l2003[ ). The flux was highly variable (by a 
factor of 20), with an usual time lapse of 10 hours be- 
tween flares, and 2 to 3 days of recurrent inactivity. The 
lines and continuum were varying on a 1000 s timescale, 
used to derive the size of the emitting region to be smaller 
than 3 x 10 13 cm. These X-ray properties, signature of 
wind accretion, were reminiscent of other peculiar high 
energy so urces, such as X TE J0421+560/CICam and 
GX 301-2 (IRevnivtsevI 12001 . 

1.1. Optical and near-infrared observations of 
IGR J16318-4848 

The accurate XMM- Newton localisat ion of 

IGR J16318-4848 allowed IFilliatre fc Chatvl (|200l 
to trigger ToO photometric and spectroscopic obser- 
vations in optical and near-infrared (NIR) just after 
the detection of the source, and to discover the coun- 
terpart at NIR wavelengths. The first striking fact 
was the extreme brightness of its NIR counterpart, 
with magnitudes J= 10.33 ± 0.14, H= 8.33 ± 0.10 
and K s — 7.19 ± 0.05. The second surprising fact 
was its absorption. The optical/NIR counterpart was 
exhibiting an unusually strong intrinsic absorption in 
the optical V-band of Ay = 17.4mags, much stronger 
than the absorption along the line of sight exhibited by 
neighbouring objects (Ay = 11.4mags), but still 100 
times lower than the absorption in X-rays. This led 
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iFilliatre fc Chatvl (|2004l) to suggest that the material 
absorbing in the X-rays had to be concentrated around 
the compact object, while the material absorbing in 
optical/NIR would extend around the whole binary 
system. 

The NIR spectroscopy in the 0.95 — 2.5/xm domain re- 
vealed the third amazing fact: the high energy source 
exhibits an unusual NIR spectrum, very rich in many 
strong emission lines. The study of these lines showed 
that they originated from different media (exhibiting var- 
ious densities and temperatures), suggesting the pres- 
ence in this high energy source of a highly complex and 
stratified circumstellar environment, and also the pres- 
ence of an envelope and a wind. Only luminous post 
main sequence stars show such extreme environments, 
the companion star being most likely a sgB[e] star, there- 
fore making IGR J16318-4848 an HMXB. As in X-rays, 
the NIR characteristics are also reminiscent of the other 
peculiar high energy source XTE J0421+560/CICam, for 
which IRAS 12-100^m data suggeste d the existence of a 
substantial circum stellar dust shell (|Belloni et al.l 119991 : 
IClark et al.lll999h . 

By fitting a black body representing the companion 
star in the multi-waveleng th Spectral Energ y Distribu- 
tion (SED) of this source, IFilliatre fc Cbatvl (|200l de- 
rived the following parameters: A\ = 17.5mags, L ~ 
10 6 L>| kpc x L , T = 20250 K, M = 30M© and r/D = 

5 x 10 -10 , where L, T, M, r and D are the companion star 
luminosity, photosphere temperature, mass, radius and 
distance, respectively. These parameters imply a high lu- 
minosity, high temperature, and high mass star, therefore 
likely a supergiant, located between 1 and 6 kpc. The 
NIR photometry, spectroscopy and SED fitting all led to 
the same results. Finally, by locating these parameters 
on a Hertzsprung-Russel (or temperature-luminosity) di- 
agram, one can see that this companion star is located 
at the edge of the blue supergiant domain, indicating 
that we are facing an ex treme object even among blue 
supergiant stars (see e.g. ILamers et al.lll998l ). 

1.2. Previous mid-infrared observations of 
IGR J 163 18- 



However, the cause of its unusually strong absorption 
could not be unveiled by observations in the optical/NIR 
domain, and only mid-infrared (MIR) data would allow 
us to characterise the nature of this absorbing material, 
and determin e whether it was mad e of cold gas, dust, or 
anything else. iRahoui et al.l (|2008l ) obtained MIR photo- 
metric observations (8- 19/xm) with VISIR on VLT/UT3 
of IGR J16318-4848. They built the optical to MIR 
SED by putting together photometric data from ESO 
NTT/Sofl, VISIR and Spitzer (GLIMPSE survey). Since 
data could not be fitted with only the spectral stellar type 
of a hot luminous sgB[e] star, they added to a model of 
the companion star (taking usual parameters of a sgB[c]) 
a simple spherical dust component. They then derived 
the following parameters: companion star temperature 
= 22000 K and radius R+ = 20 .4 R Q = 15 x 10 6 km, 
and dust component temperature Td us t = 1000 K and 
radius i? d ust = 10 R* = 150 x 10 6 km. The de- 
rived absorption and distance were Ay = 17mags and 
D = 1.6 kpc respe ctively, with a fit x 2 /dof of 6.6/6 
(|Rahoui et al.l 120081) . The most important result was 



therefore that they needed an extra component in order 
to satisfyingly fit the data, whose extension suggested 
that it was enshrouding the whole binary system, as a 
cocoon of dust would likely do. 

Archival MIR observations taken from the Spitzer 
GLIMPSE survey confirmed the existence of a long- 
wavelength ( A > 4/xm) excess in this source 

(Kapl an et al.l I20061) . These authors fit the excess with 
a black body component of temperature ~ 500 — 1500 K 
and radius ~ 10 i?*, consistent with warm circumstellar 
dust e mission (similarly to the XTE J0421+560/CI Cam 
system. rBelloni et al.lll999T ). and they suggested that this 
dust could be the cause of the X-ray absorption. How- 
ever their data could not constrain the properties of the 
dust (total mass, temperature and density distribution). 
Su bsequent Spi t zer-IR S MIR spectroscopic observations 
by IMoon et al.l (|2007D confirmed the presence of a hot 
(T>700K) circumstellar dust, and also suggested the 
presence of a warm (T^190K) dust component, which 
appeared necessary to get satisfying fits. 

Before speculating more on the nature and character- 
istics of this highly absorbed supergiant X-ray binary, we 
first need to better characterise the nature and origin of 
this surrounding dust/cold gas component. This is why 
we performed MIR spectroscopic observations of this ob- 
ject with VISIR on VLT, currently the most adequate 
instrument available to perform sensitive MIR spectro- 
scopic observations in the 7-14 /an domain, with two pri- 
mary goals: 

• to detect emission lines, in order to determine the 
composition, density and temperature of the ab- 
sorbing material. 

• to perform accurate broadband NIR-MIR SED fit- 
ting of the stellar and dust component continuum 
emission of the enshrouded binary system, in order 
to constrain the temperature, geometry, extension 
around the system, the nature and composition of 
the dust component and absorbing material. 

2. OBSERVATIONS AND DATA REDUCTION 

MIR photometric and spectroscopic observations of 
IGR J16318- 4848 were carried ou t on 2007 July 12-13 
using VISIR (jLagage et al.ll200l . the ESO/VLT MIR 
instrument, composed of an imager and a long slit- 
spectrometer covering several filters in N and Q bands 
and mounted on Unit 3 ("Melipal") of the VLT. The 
standard "chopping and nodding" MIR observational 
technique was used to suppress the background domi- 
nating at these wavelengths. Secondary mirror-chopping 
was performed in the North-South direction with an am- 
plitude of 16" at a frequency of 0.25 Hz. The nodding 
technique, needed to compensate for chopping residuals, 
was chosen as parallel to the chopping and applied using 
primary mirror offsets of 16". Because of the high ther- 
mal MIR background for ground-based observations, the 
detector integration time was set to 16 ms. 

Raw data were reduced using the dedicated IDL reduc- 
tion packageQ- The elementary images are co-added in 
real-time to obtain chopping-corrected data, then the dif- 
ferent nodding positions were combined to form the final 

3 http: / /www.eso.org/sci/facilities / paranal/instruments / visir / doc/ VISIR-dat 
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image. The VISIR detector is affected by stripes ran- 
domly triggered by some abnormal high-gain pixels. A 
dedicated destriping method was developed to suppress 
them 2 . 

2.1. VISIR Photometry 

VISIR broadband photometry was performed in six 
broad and seven narrow filters, whose characteristics are 
given in Table [TJ Each time we used the small field of 
view in all bands (19'.'2 x 19'.'2 and O'.'075/pixel plate 
scale). All the observation s were bracketed with ob- 
servations of standard stars ()Cohen et alj|1999h for flux 
calibration and PSF determination. We measured the 
instrumental fluxes and corresponding uncertainties us- 
ing aperture photometry. They were then converted to 
physical fluxes using the conversion factors derived from 
standard stars observations. Results on the photometry 
are reported in Table [1] Finally, according to the ESO 
manuaQ, systematic uncertainties of VISIR photometric 
data are of the order of 10% in PAH1 and PAH2, and 
20% in Q2. 

2.2. VISIR Spectroscopy 

We also performed VISIR low-resolution spectroscopy 

( R = ~ 35 °) of IGR J16318-4848, from 7.1 to 
13.46/im split in five filters of the N-band, centered at 
8.1, 8.8, 9.8, 11.4 and 12.4/mi. Exposure times were equal 
to 1 hour for the first filter and 2 hours for all the oth- 
ers, and we reached S/N from about 17 to 56. During 
the same night and in the same con ditions, we also ob - 
served the standard star HD 149447 (jCohen et al.lll999D . 
a K6III giant. Extraction of the spectra, wavelength 
calibration, telluric lines correction, and flux calibration 
were then carried out using the IDL reduction package 
mentioned above. In addition, in order to accurately flux 
calibrate the spectra, they have been scaled to match the 
VISIR PAH2 photometric flux. The VISIR MIR broad- 
band spectrum is shown in Figure [TJ zooms on different 
parts of the spectrum are shown in Figures [2j [3] and [4] 
respectively. The detected lines are indicated in the Fig- 
ures, and reported in Tabled 

2.3. Archival data 

With the aim of SED fitting, we completed our set of 
VISIR data with already pub lished (1) ESO NTT/SofI 
NIR (iFilliatre fc Chatvl l200l and (2) Spitzer /IBS MIR 
(|Moon et al.l 120071 ) low-resolution spectra. We have re- 
analysed and recalibrated both spectra: 

(1) The IGR J16318-4848 and HiP 80456 (F5V) tel- 
luric standard Sofl spectra, covering the spectral range 
0.9 to 2.5/im, were reduced with the IRAF suite by per- 
forming crosstalk correction, flatfielding, sky subtraction, 
and bad pixel correction. The spectra were then ex- 
tracted, wavelength calibrated with a xenon arc taken 
with the same setup, and finally combined. The telluric 
features were corrected with the telluric task. The re- 
sult was then multiplied by an F5V synthetic spectrum 
downloaded from the ESO websiteEL scaled to the 2MASS 
magnitudes of HiP 80456 in the H and K filters. The final 



IGR J16318-4848 Sofl spectrum is therefore flux cali- 
brated with about 5% uncertainties on the continuum 
level. 

(2) We used the SL2 (5.20 - 7.70/zm), SL1 
(7.40 - 14.50/zm), LL2 (14.00 - 21.30/xm), and LL1 
(19.50 — 38.00/im) low-resolution Spitzer spectra of 
IGR J16318-4848. Basic Calibration Data (BCD) were 
reduced following the standard procedure given in the 
IRS Data Handbook The basic steps were bad pixel 
correction with IRSCLEAN vl.9, sky subtraction, as well 
as extraction and calibration (wavelength and flux) of the 
spectra — with the Spectroscopic Modeling Analysis and 
Reduction Tool software SMART v8.1.2 — for each nodding 
position. Spectra were then nod-averaged to improve the 
S/N ratio. 

The broadband NIR to MIR ESO NTT/SofI, VISIR 
and Spitzer spectra are shown in Figure [SJ 

2.4. Absorption on the line-of-sight and dereddening 

The VISIR spectrum of IGR J16318-4848 shows a 
strong absorption feature, around 9.7/im, most likely 
due to silicate dust in the interstellar medium (ISM) 
along the line-of-sight. As described in Section [TJ sev- 
eral consistent measurements of the visible extinction Ay 
are found in the literature, most of them derived from 
SEP fitting, and all give values between 1 7 and 19mags 
(IFilliatre fc Chatvl[200l iMoon et alj|20"07t iRahoui et al.l 
2008) . The silicate absorption feature due to diffuse ISM 
is strongly correlated to the optical exti nction follow - 
ing the relation Ay = (18.5 ± 2) X r 9 . 7 ([Draind 12001 . 
where T9.7 is the optical depth of the silicate absorption 
at 9.70/im. It is therefore p ossible to derive it from the 
VISIR spectrum. Following IChiar et all (|2007l) . we took 
T9.7 = -ln(F 9 . 7 /F coatiimum ), where F 9 . 7 and F cont i„ UU m 
are the source and continuum fluxes at 9.70/zm respec- 
tively. The continuum was fitted using both ranges 
5.20 - 7.00^m and 13.00 - 14.50^m (in order to ex- 
clude the contribution of the silicate absorption feature) 
with a second order polynomial. We point out that 
this method is strongly uncertain, especially concern- 
ing the continuum fitting, and the result should there- 
fore be considered with caution. We nevertheless derive 
Ay = 18.3 ± 0.4mags, in agreement with the previous 
determinations. 

We used this value t o deredden the spectra with the 
extinction laws given in iChiar fc Tiele ns ( 2006). In their 
paper, these authors derived the A\/Ak ratio rather than 
the usual A\/Ay one. To express the extinction ratio in 
a standard way, we assigned to Ak the v alue derived 
from the law for the diffuse ISM (Ry = 3.1. lFitzpatrickl 
119991) : A K = 0.111 x Ay Moreover, these authors find 
that beyond 8/um, the MIR extinction and silicate ab- 
sorption features are best described by two methods, us- 
ing: 1) the diffuse ISM (which stops at 27/xm), and 2) 
the Galactic center. We used both cases to deredden 
our spectra of IGR J16318— 4848, however the Galactic 
center case clearly overestimates the silicate absorption 
depth. We therefore adopt the corrections based on the 
extinction due to the diffuse ISM. 

3. RESULTS 



4 www.eso.org/ sci / facilities / par anal/instruments / visir/doc 

5 http: / /www. eso.org/sci/facilitics/paranal/instruments/isaac/tools/lib/lrlrfibp:^Bsh. spitzer. caltech.edu/irs/irsinstrumenthandbook/IRS_Instrument_F 
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3.1. Photometry 

We first compare in Table Q] IGR J16318-4848 MIR 
photom etric fluxes with ou r two previous VISIR obser- 
vations ()Rahoui et al.l l2008). to see whether there is vari- 
ability in the MIR emission of this source. By examining 
these values, it appears that there exists a small variation 
of the fluxes: 

• the PAH1 flux has increased by ~ 6 — 7a between 

2005 and 2006, then remained constant taking into 
account the error bars between 2006 and 2007. 

• the PAH2 flux has decreased by 1.5a between 2005 
and 2006, then further decreased by ~ 8a between 

2006 and 2007. 

• the Q2 flux remained constant taking into account 
the error bars between 2005 and 2006, then in- 
creased by ~ 2a between 2006 and 2007. 

There are therefore 2 significant (> 6a) variations of 
the flux: first the increase of the PAH1 flux between 
2005 and 2006, then the decrease of the PAH2 flux be- 
tween 2006 and 2007. All the other variations are not 
significant (< 2a). Con side ring the systematic uncer- 
tainties given in Section 12.11 we consider in this paper 
that these variations are within the error range and that 
the source remains at a stable flux. We accordingly use 
this photometry to better calibrate the spectra, as de- 
scribed in Section 12.21 and shown in Figure El in order to 
accurately compute fits of IGR J16318-4848 SED. Fu- 
ture additional photometric observations should allow us 
to better constrain the potential flux variability on a long 
timescale. 

3.2. Spectroscopy 

With the spectroscopic observations, we first confirm 
the extremely rich and complic ated MIR environm ent of 
the source already reported bv lMoon et al.1 (|2007| ). how- 
ever we detect many more emission lines in our MIR 
spectra obtained with the VISIR instrument. In these 
spectra, presented in Figures [TJ O El HI and also [5j we 
clearly see the following lines, reported in detail in Table 

m 

• the broad and prominent silicate absorption fea- 
ture around 9.7/im, due to dust c omponent in the 
interstellar medium (|Drainell2003f ); 

• numerous H I recombination lines - the most 
prominent being H 1(7-6) at 12.35/im - due to ion- 
ized stellar winds; 

• many forbidden atomic and ionic fine-structure 
lines of low ionization potential in emission 
Ne II], [Ne VI], [Na IV], [S IV], [Ar III], [Ar V], 
Ca V], [Ni II], and tentatively [Co II]-, caused by 

extended ionized regions of low density; 

• many higher ionization He II emission lines; 

• and various polycyclic aromatic hydrocarbon 
(PAH) emission lines, originating from photodis- 
sociated regions. 



Some lines, especially H and He emission lines, exhibit 
P Cygni profile. The presence of many forbidden lines 
of low excitation metals such as [Ni II] indicates that the 
circumstellar material is geometric ally extended, with a 
large amount of low density gas (jLamers "et~atl[T998h . 
Mos t of the lines that w e report exhibit the same fluxes 
as in lMoon et all (|2007l ). apart from H 1(16-9) and [Ni II] 
(fainter by a factor 2) and the [Ne II] 12.8jum li n e (ten 
times fainter). In the spectrum of iMoon" eTaTl (120071) 
there is a line at the position of [Ne II] which is much 
stronger than that in our spectrum, however since their 
spectral resolution is lower than ours, it is difficult to 
compare fluxes of these lines. On the other hand, X-ray 
emission coming from the compact object might pho- 
toionize and increase the [Ne II] 12.8/im line flux, and 
not higher excitation states. 

We can use various lines of the MIR domain as indi- 
cators of the physical medium state of the source, sensi- 
tive to t he temperatur e of th e environment, based on the 
work of iFurness et al.l (120101 ) (to our knowledge there is 
no other line diagnostics that we could use in this MIR 
domain). Since for most of the lines there is no significant 
variation, we also use for these diagno stics two additional 
lines coming from IMoon et al.l ((2007): [Nc III] 15.56/xm 
and [S III] 18.71/zm. We therefore find that: 

, , [Aril I] 8.99/mi x _ 
l0 9\ [NeII\ 12.81/inV — 

-0.029, indicates an effective 
38000- 40000 K. 



• The intensity ratio 

<2.02e-21 



log{\ 



21 > 



^2.16e-21 

temperature T e g 



• The intensity 
perature T e ff « 



ratio lnn( [NeIII] 15 - 56 ^ n V - 

ratio, Log{ [NeII] 12 .sx/dn > ~ 
0.08, indicates an effective tem- 
38000 - 40000 K. This intensity 



ratio correlates well with log( [^e//( 12 si^m )' tak- 
ing into account the scattering. If inst ead we take 
the va lue of [Ne II] 12.81/im given in IMoon et al.l 
(|2007f ). and therefore obtained at the same time 
than [Ne III] 15.56/im, we derive a value of ~ 
35000- 36000 K. 

• The intensity ratio, log( ) = 

^0( 2 16e-2i ) = — 0.29, indicates an effective 
temperature T cff ~ 37000 - 38000 K. This ratio 
agrees well with log( TC^ryrr Tjwfifir )' taking into 
account the scattering. 

• The intensity ratio, I og(j§jjjp^^) = 
log( ifUfEf? ) — —0.84, is close to the correla- 
tion with log( p^r^r w sijtm )' taking into account 
the scattering and uncertainties. 

• the log(r) < S- Ne>) versus log( fo 6 //^ l^si^ ) 
corresponds well to an effective temperature T c g ~ 
38000 K with l ogU = — 1, adopting the same no- 
tations than in IFurness et al.l ([201Q;), in particular 

n <r ^-Np ( \- SIV \ 10-52/xm ]/( [NeIII] 15.56 M ni x 

// ^ a iVd ^— V[s/JJ] 18.71/mi" V [Nell] 12.81pm /' 

As discussed later in Section H this ~ 37000- 40000 K 
temperature -higher than a ~ 20000 K photosphere of a 
BI star- corresponds to a heated zone surrounding the 
stellar photosphere. 
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3.3. Spectral energy distribution fitting 

The presence of PAH and fine-structure emission lines 
such as [Ne II] or [Ar III] proves the existence of a pho- 
toionized dust component intri nsic t o IGR J16318— 4848 . 
Moreover. iKaplan et al.l ([20061 ) and lRahoui et all (2008) 
detect ed strong NIR and MIR excesses, usual in sgB[e] 
stars (jLamers et al.1 119981 ) , that they fitted with a ~ 
1000 K spherical black body, likely corresponding to 
warm circumstellar dust, with an extension of ~ 10 i?*. 
Nevertheless, this modeling is too simplistic, for the 
two following rea sons. First because, as shown by 
iMoon et al.l j2007) , by fitting the spectral continuum 
from 5 to 35/im instead of photometric fluxes, a cold ~ 
190 K dust component is needed in addition to the warm 
one. Second because, if the IGR J 163 18— 4848 compan- 
ion sta r is indeed a sgB [e] as shown by iFilliatre fc Chatyl 
(2004), the dust responsible for the NIR/MIR excess 
should thereby be distributed in an equatorial disk-like 
str ucture, instead of pre senting a spherical geometry (see 
e.g. iLamers et al.lll998f ). 

Recently, several improvements have been achieved 
in the understanding of the emission of intermediate- 
mass stars such as Herbig Ae/Be stars (HAEBE). Be- 
sides the star, these systems are known to exhibit an 
equatorial disk composed of (see Figure [6]): i. an in- 
ner dust-free cavity filled by optically thin gas, ii. a 
hot and optically thick irradiated rim made of puffed-up 
hot dust at the sublimation temperature - responsible 
for the NIR excess-, and Hi. a surrounding component 
of warm dust shell responsible for the MIR ex cess (see 
e.g. iDullemond et al.ll200ll iMonnier et~aTll2005l and ref- 
erences therein). 

This scheme was simplified bv lLachaume et all ((2007), 
by assuming that the disk was completely flat, with a 
ring- like rim of constant temperature T r j m , located at an 
inner radius R- ln and of width -ff r j m (see Figure [6]) . The 
remaining parts of this disc are located at an inner ra- 
dius i?i n + -ffrim/2, modelled with a radial temperature 
profile T (R) as given below in Equation [2] In the case 
of HAEBE, the rim is assumed to be located at the dust 
sublimation radius, with a temperature of 1500—2000 K, 
in order to put a constraint on R m . In the case of sgB [e] 
stars the presence of strong O I emission lines detected in 
their NIR spectra, as it i s the case for IGR J16318-4848 
(jFilliatre fe Chatvll2004h . rather points to wards a tem- 
perat ure between ~ 5000 and 10000 K (jKraus et al.1 
2D07). Therefore, since there is no constraint on the disk 
inner radius R m value, it will be a free parameter, to 
make the fit more accurate. We adopt the emission model 
summarised in Equation[2]to fit IGR J16318-4848 SED: 

F v = (^-) x B{y, T.) + 2 Hti ™f n cos (i) B{v, T^Q) 

+2 CO£j>) ( R ° M RB ( v T ( R))dR 
D* JR in 

with T(R) = T in (R/R in )t2) 

where B is the black body function at the frequency v, 
T*, i?* and the star temperature, radius and distance 
respectively, and i, Ti„, i? ut and q the disk inclination 
angle, inner region temperature, outer radius and a di- 
mensionless "temperature index" parameter respectively. 



T and R are the generic temperature and radius respec- 
tively. 

We can re- write Eq. [5] using the new variables r* = jf-, 
a generic radius r = R ^^}^ , r in = R in ^ c ° s ^ , r ou t = 

a/ cos(i) , TT -J cos(i) , , . 

Rout^^i — , h lim = i/rim^-o; — , and we obtain: 

F v = r* 2 B{v,T*) + 2h llm r in B(v,T I Q) 
+2 / rB{v,T{r))dr 

with T(r) = T in (r/r in )- q (4) 

During the fitting process, we fixed T„ = 20000 K 
(consiste nt with the temperature of a sgB[e] as de- 
rived by IFilliatre fc Chatvl [20041 . iRahoui et all [20081) 
and considered two cases for the temperature index 
parameter: q = .5 for an irradiated/flaring disk 
(IChiang fc Goldreichl[T997h. an d q = 0.75 for a viscous 
disk (jShakura fc Sunvaevlll973h . 

We first fit the data using the flaring disk model, how- 
ever in each case but one (reported in Table [3]) , Trim is 
found to be lower than the inner dusty disk temperature 
Tin, and we can not obtain any constraint on a param- 
eter set corresponding to high T r - lm . This result does 
not appear physical since it is not consistent with sgB[e] 
characteristics, and it likely means that the flaring disk 
model is not adapted for a sgB[e] star. 

We then fit the data using the viscous disk model, and 
the best fit is obtained for an irradiated rim of tempera- 
ture T r i m = 3786 K, along with a hot dusty viscous disk 
component at T in = 767 K. The derived best-fit param- 
eters are listed in Table OH and the broadband NIR to 
MIR resulting fit on ESO NTT/SofI, VISIR and Spitzer 
data is reported in Figure CO 

However, since this value seems quite low for a sgB[e] 
star exhibiting O I emission lines (see above), we tried 
various additional fits by freezing reasonable T r i m going 
from 4000 to 6000 K. For T rim higher than 5963 K, the 
stellar contribution is not necessary anymore, which is 
very likely unreal, since we clearly see its contribution 
in the NIR domain. We also point out that all the fits 
with T r i m between 3786 K and 5500 K pass through the 
I band data point (taking into account the error bar), 
which is not the case for T r - lm higher than this. This is 
another argument in favour of T r i m in this range. For 
these cases, the most satisfying fit seems to be reached 
for an irradiated rim of temperature T Y - lm ~ 5500 K, along 
with a hot dusty viscous disk component at T- ln = 895 K. 
The derived parameters are listed in Table El and the 
broadband NIR to MIR resulting fit on ESO NTT/SofI, 
VISIR and Spitzer data is represented in Figure [5] 

Finally, we also compare fits using i. Sofl-VISIR- 
Spitzer data and ii. Sofl-Spitzer data (i.e. without 
VISIR data). First, by taking components of spheri- 
cal geometry only, both fits exhibit similar goodness, 
with SofI-VISIR-5piteer data fits leading to a slightly 
smaller reduced x 2 of 1.29, compared to 1.32 for Sofl- 
Spitzer data, and best fit dust temperatures of 2668- 
3108 and 379-412 K for both spherical components re- 
spectively (see Tables [3] and [4] at the "2-sphere compo- 
nent" line). In addition to a high value of the reduced 
X 2 , it is the temperature of the first spherical component, 
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too high to be consistent with any kind of dust, which 
tends to exclude the spherical alternative. Second, and 
this is the most interesting result, by taking components 
of disk-like geometry, the fit is significantly better with 
Sofl-VISIR- Spitzer data (reduced \ 2 °f 1-19 compared 
to 1.30 for SoR-Spitzer data, with best fit temperatures 
of 3786 and 767 K for both components respectively), 
probably due to the fact that more data points are avail- 
able in this crucial wavelength range. Therefore, while 
Spitzer-oiily data do not allow us to distinguish between 
the two geometries, adding VISIR data lead to a quan- 
titative fit clearly in favour of the disk-like geometry. 
By excluding the spherical geometry for the irradiated 
component, VISIR data are thus essential for a better 
assessment of the dust component geometry surrounding 
IGR J16318— 4848. 

4. DISCUSSION AND PERSPECTIVES 

The multi-wavelength observations that we have per- 
formed from optical to MIR on the highly absorbed 
HMXB IGR J16318-4848 hosting a supergiant B[e] star 
allow us to strengthen the presence of absorbing material 
(dust and cold gas) enshrouding the whole binary sys- 
tem and better characterise its properties. With these 
new MIR spectroscopic observations and NIR-MIR ESO 
NTT/Sofl, VISIR and Spitzer spectral fitting we confirm 
the presence of two-temperature components, with an ir- 
radiated rim of temperature Ti; m = 3786 K, surrounded 
by a hot dusty viscous disk component at T- ln = 767 K, 
corresponding to th e best fit paramete r s (see Table |3j) . 

As sug gested by IFilliatre fc Chatvl (120041 ) and then 
shown bv Kaplan et al l (12006ft. iM oon et al.l (|2007l) and 
iRahoui et al.l (|2008D . our MIR spectra and SED fitting 
confirm the existence of a strong NIR/MIR excess, due to 
the presence of a hot circumstellar dust component, likely 
formed in the dense outflows coming from the early-type 
sgB[e] companion star of IGR J16318— 4848. This com- 
ponent has a temperature of Td us t ~ 767 — 923 K, with 
a maximum extension of r out = 3.55 au/kpc (assuming 
cos(i) = 1) for R+l D+ betwe en 5.80 - 11.71 iW kpc. 
Assuming a distance of 1.6 kpc (Ra houi et alJl2"008ft . i?* 
would be in the range 9.3 — 18.7 Rq, and the exten- 
sion of the dust equal to r out = 5.6 au, corresponding 
to 64.3 — 129.8 R+. This value is larger than the one de- 
rived with the spherical shell model, which is consistent 
with the fact that for the same amount of dust, it will 
have to be spread on a larger distance if it has a disk 
rather than shell geometry. In addition, a high value 
is consistent with a low inclination, a s suggested by the 
Br7 l ine profile in the NIR spectrum (jFilliatre fc Chatvl 
12004ft . 

On the other hand, in none of our NIR-MIR fits do 
we require a warm ~ 190 K dust component that was 
claimed by IMoon et all (|2007| ) . We suggest that this ad- 
ditional component was necessary in their fits, firstly be- 
cause they only took into account the Spitzer spectrum 
without including the NIR spectrum, and secondly be- 
cause they chose a spherical geometry. By including all 
available spectra from NIR to MIR, and applying a disk- 
geometry to our components, our fits are much closer to 
the physical nature of the environment surrounding these 
objects, and therefore more accurate and physical. 

We will now derive the mass of the dust which is 
hosted in the viscous disk surrounding the sgB[e] star. 



By extrapolating our best- fit spectrum (see Figure [7]), we 
should have F v = 0.155 Jy at 60/im. At this wavelength, 
taking an initial MRN distr ibution with thin i ce mantles 
(Mathis Rumpl Nordsieck; iMathis et ai]|1977j), the ISM 
dust opacity is «v(60/zm) ~ 87 cm 2 /g (|Drainell2003ft . We 
compute the black body function for T; n = 767 K, corre- 
sponding to the dust inner temperature given by the best 
fit (see Table EJ): B{v, T dust ) = 1.57 x 10 15 Jy. Then, 

F 

the total dust mass is given by Mdust = ~ B { v T * n 
where F v — 2 J^° ut rB{v,T (r))dr, with r- m and r out cor- 
responding to the inner and outer radii given by the best 
fit, respectively equal to 0.74 and 3.47|^ y / (cos(i)) (see 
Table [3|). We therefore derive a dust mass in the disk of 
M dust ~ 5.43 x 10~ 9 An k p C 2 M Q , where D+ lk is the 
star distance expres sed in kpc. For a distance of 1.6 kpc 
(|Rahoui et al.ll200l . we find M dust - 1.39 x 10~ 8 M Q = 
2.76 x 10 22 kg. 

Thi s dust mass is lower than derived by IMoon et al.l 
( 2007) (3.1 x 10" 8 M Q and 4.6 x 10~ 10 M for the warm 
and hot components respectively). This is consistent 
with the fact that for a same obscuration, a shell ge- 
ometry will require more materi al than a disk geom etry. 
This also shows that, unlike in Mo on et al.l (|2007[ ). our 
fits do not require any hot component. 

There are a few possibilities to explain the low mass 
of dust composing this stellar disk, compared to e.g. 
standard giant stars surrounded by disk, where the dust 
mass in the disk is more o f the order of 10 29 g (see e.g. 
Uura et al.ll200lH Jurall20~0l : i. the disk might be smaller 
in sgB[e] stars due to evolution effect, ii. special condi- 
tions such as X-ray ionisation might destroy the dust in 
the disk, and/or Hi. the compact object, accreting from 
the stellar wind, might have an influence on the mass 
of dust comprised in the disk (since the compact object 
likely orbits inside the disk, as suggested by the high 
absorption and the extension of the dust). 

The low mass of dust could also be explained by 
tidal torques created by the star, disk and compact ob- 
j ect, and causing the disk t o be potentially truncated 
(|Okazaki &: Neg ucrucla 2001), likely at the n n distance. 
The compact object might orbit within or at the edge 
of the disk, with a small inclination, quasi-coplanar with 
the disk, causing the extreme absorption seen in X-rays, 
probably local around the compact object (see Section!!}. 
In this case the wind accretion would be fully consistent 
with the disk geometry derived by SED fitting. However 
it is difficult to conclude on this disk truncation, since 
the orbital separation is unknown for this system. 

All the line ratios reported in Section 13.21 indicate an 
effective temperature T cff ~ 35000 - 40000 K. This ef- 
fective temperature is higher than the temperature of a 
hot Bl supergiant sta r, of luminosity between 10 5 and 
10 6 L Q . For instance, IFilliatre fc Chatvl (|200l derived 
a temperature of T = 20250 K by fitting the optical-NIR 
SED. However we point out that while the NIR spectrum 
mainly emanates from the stellar photosphere, the MIR 
spectrum comes from the heated rim and the dusty vis- 
cous disk, as is shown by both the broadband and best- 
fit spectra (see Figures [5] and |7J). Therefore the location 
where the MIR lines are produced does not correspond to 
the stellar photosphere, but instead to dense and hot re- 
gions protected from the stellar radiation, and probably 
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irradiated by X-ray emission coming from the compact 
object accreting from the stellar wind. In addition, this 
X-ray irradiation could explain the variability of a few 
lines, while most lines exhibit a stable flux. 

We used in the c omputations a bove a distance of 
1.6 kpc as derived bv iRahoui et al.l (|2008| ). however this 
is the value derived with a star dominating at NIR wave- 
lengths, and it is possible to have a star located further 
away. The best fit gives R+/ D+ = 11.71 Rq /kpc. A 
distance of 1.5 — 2 kpc gives R± = 17.6 — 23.4 Rq, corre- 
sponding to a typical radius of supergiant star. However, 
considering the more physical fit with T r i m = 5500 K, the 
best-fit parameter R+/ = 5.80 Rq /kpc requires a 
distance of 3 — 4 kpc to get a stellar radius correspond- 
ing to a supergiant star (see Table [3j> - 

In the broadband spectrum we see two additional ab- 
sorption components, at around 10 and 20/xm (see Figure 
likely corresponding to silicate absorption. We point 
out that these components might be due either to an ab- 
sorption law which is not accurate enough in this MIR 
domain, or even to auto-absorption of a disk-like system 
that would be seen nearly edge-on, which in this case 
would be suggesting a high inclination system. 

These observations allowed us to relate the presence of 
the hot dust component to the nature of the early-type 
star and to stellar wind photoionization. Indeed, sgB[e] 
stars usually exhibit hot (T~1000K) circumstellar dust, 
and show forbidden line emission in the MIR (see e.g. 
lLamers et al.l 119981 ) . While the evolutionary status is 
not well known yet, they will likely evolve into extreme 
and still poorly known Luminous Blue Variables (LBV) 
or Wolf-Rayet (WR) stars. MIR spectroscopic character- 
istics of IGR J16318— 4848 appear very similar to those 
from LBVs, s ome of them exhibit ing also thermal MIR 
dust emission (jLamers et al.lfl996l ) . If the evolution from 
sgB[e] to LBV is confirmed, these highly obscured IN- 
TEGRAL sources might therefore be the progenitors of 
LBV stars, or WR stars in case of high mass loss, or even 
dusty wind hypergiant stars. If this is the case, then ob- 
scured INTEGRAL sources might represent a previously 
unknown evolutionary phase of X-ray sou rces hosting 
early -type companion stars (see discussion in lMoon et al.1 
2007). 

The question then arises: does the supergiant star 
influence its close environment, or does the interstellar 
medium influence the star, by a feedback effect? As dust 
is a strong tracer of star formation, the next step is to 
observe these sources with Herschel, since its greatest 
strength is the possibility to study the history of star for- 
mation in our Galaxy. This makes observations of such 
obscured sources with Herschel an enormously powerful 
tool, giving a new insight for studying the process of 
formation and evolution of very massive stars in molec- 
ular clouds. These observations should be valuable in 
the study of the later phases of stellar evolution, particu- 
larly circumstellar shells, mass-loss in general, and stellar 
winds, and they should also allow us to relate the proper- 
ties of these objects to their birth environment. Higher 
wavelength observations should also allow us to better 
constrain the broadband SED and to confirm T r j m and 
Ti n we derived from our observations. They would likely 
close the debate on the presence of colder (< 200 K) 
circumstellar dust component around this source, with 
an extension at larger radii from the compact object of 



more than 30 i?*, suggested by Mo on et all (|2007ft but 
not confirmed by our observations. Finally, these obser- 
vations should allow us to investigate where this material 
comes from, and the study of such an unusual circumstel- 
lar environment should in the future tell us whether it is 
due to stellar evolution or to the binarity of the system 
itself. 

5. CONCLUSION 

We have performed new VISIR photometric and spec- 
troscopic MIR observations of the highly absorbed 
HMXB IGR J16318-4848 hosting a supergiant B[c] 
star, and we have fit its NIR-MIR SED including ESO 
NTT/Sofl, VISIR and Spitzer data. The main results 
are: 

• the photometry did show two periods of significant 
variability: first an increase in PAH1 between 2005 
and 2006, then a decrease in PAH2 between 2006 
and 2007. Apart from this the source remained at 
a constant flux, within uncertainties. 

• the spectra exhibit many emission lines from vari- 
ous elements, some coming from forbidden atomic 
and ionic fine-structure lines of low ionization po- 
tential. The line diagnostics all point towards a 
source of temperature between 35000 and 40000 K, 
likely emanating from a heated zone surrounding 
the stellar photosphere. 

• the multi-wavelength optical to MIR Sofl-VISIR- 
Spitzer SED fitting allowed us i. to exclude a 
spherical geometry for the dust component, and ii. 
by using a model adapted for sgB[e] stars, to show 
the presence of two-temperature components, with 
an irradiated rim of temperature T r i m = 3786 K, 
surrounded by a hot dusty viscous disk component 
at T in = 767 K, corresponding to the best fit pa- 
rameters (see Table [3l). However, the presence of 
O I lines rather points towards irradiated rim tem- 
peratures of ~ 5500 K, for which we achieved sim- 
ilar goodness of fits. 

These results allowed us to strengthen the presence of 
absorbing material (dust and cold gas) enshrouding the 
whole binary system of IGR J16318— 4848, and better 
characterise its properties. They are therefore of prime 
importance to constrain the nature of this system in par- 
ticular. In addition, they will also allow us to better 
understand the whole population of new high energy bi- 
nary systems born with two very massive components, 
which has been recently revealed by INTEGRAL, these 
sources being the most extremely absorbed ones among 
high energy sources. These features are most likely asso- 
ciated with the mass losses of OB supergiants and may 
play an important role in their evolutions to the final 
supernova explosion and formation of compact objects. 
These systems are probably the primary progenitors of 
neutron star/black hole binary mergers, with the pos- 
sibility that they are related with short/hard gamma- 
ray bursts, and that they might be good candidates of 
gravitational wave emitters. It therefore appears that, 
because these obscured X-ray binaries represent a differ- 
ent evolutionary state of X-ray binaries, their study is of 
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prime importance, and will provide hints on the forma- 
tion and evolution of high energy binary systems, and a 
better understanding of the evolution of such rare and 
short-living HMXBs. 
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Fig. 1— Broadband ESO/VLT VISIR MIR spectrum of IGR J16318-4848 from 7.2 to 13.6 [tm. Particularly visible is the broad and 
prominent silicate absorption feature around 9.7/im, due to dust component in the interstellar medium. 
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Fig. 2.— Enlarged ESO/VLT VISIR MIR spectrum of IGR J16318-4848 from 7.2 to 9.3 /»m with identified lines indicated. 
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Fig. 3.— Enlarged ESO/VLT VISIR MIR spectrum of IGR J16318-4848 from 9.3 to 11.4 ^m with identified lines indicated. 
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Fig. 5.— Broadband NIR to MIR ESO NTT/Sofl, VISIR and Spitzer reddened spectrum of IGR J16318-4848, from 0.9 to 35/xm. The 
"corrected photometry" corresponds to VISIR + Spitzer photometric data, taking into account the transmission factor. 
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Fig. 6. — Representation of the sgB[e] sche me, with the 3 components seen from polar view (this diagram is just for illustration purpose, 
and is not to scale, see details in Section l3.3|l . 
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Fig. 7.— Broadband NIR to MIR fit on ESO NTT/SofI, VISIR and Spitzer dereddened spectrum of IGR J16318-4848. SED with a 
stellar spectral type corresponding to a sgB[e] with T, = 20000 K, T Ilm = 3786 K and d ust temperature of ~ 767 K. These parameters 
correspond to the best fit, however T r i m seems very low for a sgB[e] star (see text, Section 1 3 - 3 P> - 
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Fig. 8.— Broadband NIR to MIR fit on ESO NTT/SofI, VISIR and Spitzer dereddened spectrum of IGR J16318-4848. SED with a 
stellar spectral type corresponding to a sgB[e] with T* = 20000 K, T r i m = 5500 K and dust temperature of ~ 895 K. This is not the best 
fit, however it seems physically more real to describe a sgB[e] star (see text, section [3. 3 P> - 
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TABLE 1 

VISIR Photometry of IGR J16318-4848 a . 



Filters X(pm) AA(yum) Flux 2005 (mjy) Flux 2006 (mjy) Flux 2007 (mjy) 



PAH1 


8.59 


0.42 


409.2 ± 2.4 


426.2 ±3.0 


420.98 ± 3.07 


Arlll 


8.99 


0.14 






281.25 ±4.92 


SIV.l 


9.82 


0.18 






210.93 ± 14.79 


SIV 


10.49 


0.16 






249.63 ±4.12 


SIV.2 


10.77 


0.19 






254.59 ±4.22 


PAH2 


11.25 


0.59 


322.4 ± 3.3 


317.4 ±3.4 


290.96 ±3.01 


SIC 


11.85 


2.34 






361.33 ±3.16 


PAH2.2 


11.88 


0.37 






321.02 ±3.39 


Nell.l 


12.27 


0.18 






372.69 ± 5.75 


Nell 


12.81 


0.21 






376.74 ±8.31 


NeII.2 


13.04 


0.22 






365.98 ± 7.38 


Ql 


17.65 


0.83 






232.24 ± 14.40 


Q2 


18.72 


0.88 


172.1 ± 14.9 


180.7 ± 15.3 


208.13 ± 11.27 



a Photometry taken on 2005 June 21, 2006 June 30 dRahoui et alj 120081) and 2007 July 12 
(this paper) . 
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TABLE 2 

Detected lines in the IGR J 16318-4848 VISIR MIR spectrum 



Identification A (/im) A flt (/urn) Flux (10- 21 W/cm 2 ) SNR flt EW (/urn) FWHM (10- 2 /zm) 



HI (6-5) 


7.4599 


7 4638+2 4361 X 1 0" 

1 m HJtJO _|_ Z . TTtJ VJ 1 /X 1U 


-4 


824 08+39 37 


30.17 


-0 51 16.3 


1 5894+2 4401x10" 

1 . juci _i z . uu 1 /\ 1 vj 


-2 


[NcVI] 


7.6524 


7 6630+1 0785x10" 

1 . YJ\J*-JVJ _l L.U 1 OcJ /X 1U 


■3 


148 91+32 60 


6.86 


-0 06377 

VJ 1 VJVJ tJ 1 1 


1 5871+1 1351 XlO" 

1 . tJVJ 1 1 1 L. IrJtJl /X 1VJ 


-2 


PAH 


7.7000 


7 7080+5 6055 XlO" 

1 . i VJOVJ _1 7 . V J V ' - ' 7 /X 1U 


-4 


27 67+2 74 

^ 1 >U 1 _l_ ^. 1 1 


15.64 


-0 00924 

vj . vjvj c z/i 


1 8649+6 3628x10" 

1 . O VJ iJ 1 VJ . tJ VJ Z t J» /X 1 VJ 


-2 


HT Clfi-81 

111 1 1U O J 


7.7804 


7 7511+9 0737x10" 

1 i I Jl J — 1_ C . W f » J 1 /X 1 VJ 


-4 


44 00+6 00 


9.20 


-0.01771 


2 0798+9 0719x10" 

z . vj 1 c<j 1 \J . vj 1 u a iw 


-2 




7.9016 


7 8706+1 8145x10" 

1 .(3 1 VJVJ _l L ,011'J /x 1U 


3 


74 64+1 7 1 fi 

1 1-U1J L f . 1 VJ 


4.95 


-0 03614 

VJ . UiJVjH 


2 c;016+1 8166x10" 

Z . • J VJ 1 VJ _l 1.01 V J VJ /x 1VJ 


-2 


[rtl V J . 


7.9016 


7 9274+3 8353x10" 


-4 


43 67+3 72 


17.95 


-0 01898 

VJ . W _L O t J O 


1 3931+3 8038x10" 

1 i'jj jj — 1 .J . vj » J v ' /x 1 vj 


-2 


HoTT ("21-141 /Ho 

1 1 V>1 1 \ Z 1 11 1 / 11^ 


8 0385 


8 0309+9 4853x10" 


-4 


16 26+3 30 


7.53 


-0 00764 

vj • vjvj j v ' 1 


1 4821+9 7625x10" 

l.lUil — 1_ C . 1 UZ vj /X 1 VJ 


-2 


HpTT ("24-1 5ln- 

1 1V_>1 1 1 Zrfc 1 vj T U 


8.4128 


8 ^882+1 0484x10" 


■3 


k qi -i-o 90 


7.73 


+0 00330 

1 VJ . VJ VJ- J tJVJ 


1 0228+I 0491x10" 

1 . >jZZO_l L . VJ'-tiJ 1 A 1U 


-2 


HeTT ("24-15ln+ 

11V>11 l Zrt ly J |_/ | 


8.4128 


8 4187+5 7743x10" 

O . H O 1 _l_ ".J . ■ t 'J s\ _L W 


-4 


3 63+1 05 


9.38 


-0 00248 

VJ . VJVJZrtVJ 


95614+9 3446x10 

u.tj!Jum_Ct(jnvj a ivj 


-2 


HT (25- Q 1r>- 

ill 1 ZvJ ^ J J-J 


8.4849 


8 4772+1 1904x10" 

O . 1 1 1 Z _l L. UUT A 1 VJ 


■3 


7 03+1 42 


6.66 


+0 00499 


1 8254+1 1948x10" 

1 . uzuu l . 1 iU /x ivj 


-2 


HT f25-9ln4- 

111 1 Z/O J /pT 


8.4849 


8 5066+7 8518x10" 


-4 


k e;o+2 52 


6.48 


-0 00466 

VJ . VJVJUJVJ 


95303+1 4639x10 

VJ . tJ'-JiJVJcJ—l LWOJA IVJ 


-2 


PAH 


8 6000 


8.6216+1.1826x10" 


3 


7.70+1.34 


7.10 


-0 00534 

vj . v ' v 7 - ' ■ J 1 


2.1248+1.1826x10" 


-2 


HI (14-8) 


8 6645 

O . VJVJIvJ 


8 6844+1 8506x10" 

O . VJOH—I I . ( ; : / v » V ' /\ 1 V ' 


■3 


4K-I-1 82 


6.95 


-0 0071 5 

VJ . VJVJ 1 ItJ 


•3 18144-2 2038x10" 

■J . 1 H 1 Z . Z VJ » J O -^x IVJ 


-2 


HT ?23-Ql 

111 \ ZtJ £7 1 


8.7206 


8 7iqR-|-fi ^877x10" 

U. 1 1 ■ " J 1 1 J . <JO 1 1 /x 1 VJ 


-4 


09 09-1-1 07 

Zj Zj . _l L • i7 ■ 


13.84 


-0 02097 

VJ . V J Z V J iy 1 


2 41 q^+fi 706"^ x 1 0" 

Z . 1 1 C7 !j _l_ VJ . 1 VJVJ ij /x IVJ 


-2 


HT ("10-7ln- 

111 \ IVJ 1 J \J 


8.7601 


8 7fi1f|-(-7 19S1 xlO" 

O . 1 U1U _l_ 1 .UOl A 1U 


-4 


1 n 86+1 29 

1 VJ . O VJ _l L . Zj rj 


11.21 


+0 01 003 

I VJ . VJ 1 VJVJ'J 


1 8906+7 2333x10" 

1 . OiJVJVJ _l_ 1 .ZiJtJiJAlVJ 


-2 


HT f10-71n-l- 


8.7601 


8 701R-I-2 ^839x10" 

O. 1 tJ 1VJ_I_Z . * J O KJ U /x 1 VJ 


-4 


1 8 08+1 1 2 

IO . £70_l L.lZi 


25.52 


-0 01 857 

VJ . UlOw 1 


1 QQ78+2 471 5x10" 

1 . > 1 O _l_ ii.1l lij A 1U 


-2 


HT C22-Q1 

111 I Z Z tj J 


8.8697 


8 8726+3 2955x10" 

O . O 1 Z VJ _l_ -J . Z vJO tJ /x IVJ 


3 


•3 25+1 04 

tJ . iiO_l L . VJ'i 


3.28 


-0 00320 

VJ . VJVJtJZVJ 


•3 Qqsi -i-Q ^891 x 1 0" 

•J . > J iy ' > J — L'J . Jui? 1 A ivj 


-2 


HpTT ("23-1 51 

llV^ll I^O 1 J| 


8 9209 

O . CZUy 


8 qi Qq4- rr 8644 x 1 0" 

O . ClUf _l_ ;J . OVJH A 1U 


-4 


1 88+0 26 


11.25 


-0 001 89 

VJ . VJVJ 103 


1 2705+5 6207x10" 

1 . Z 1 VJtJ_l_!j . VJZVJ 1 A 1U 


-2 


[Arlll] 


8.9914 


o.yooon.oi loAiu 


-3 


9 09-1-0 


7.43 


001 Q1 
u . uu 1 y 1 


^ 40^17-1-9 141^v10~ 

O . 1Uu 1 . Irt IO A 1U 


-2 


fNaTVl 

1 N dl V 


Q ozn oo 


y.UOOIIU.UUU 1 ! A ±KJ 


-4 




10.62 


0090fi 

-U . UUzUU 


1 4c:qi _!_(? fi^41v10~ 
1 . uoi _r_u . uoii a iu 


-2 


HT ("20-91 

111 1 ZVJ iJ 1 


vj . Z VJ VJ !j 


. Zj'-t'J VJ _l_ 1 .Oi/IJAIU 


-4 


2 ^2+0 *3(S 

Zj . O Zj _l_ V J . • J V J 


9.44 


VJ ■ VJVJcJ' J C7 


1 «qqc:_|_o 0Q2Qy10" 


-2 


HT ?1 3-81 

111 1 1J O J 


9 3920 

v? . cJiJZVJ 


n •3701-1-1 0614x10" 

i/.O 1 £71 _l L . \J\J 11 A 1U 


-3 


r> 8^ + 80 


7.09 


-0 00723 

VJ . VJVJ 1 ZtJ 


1 koq7_|_i 0270x10" 

1 . uJJ 1 _l L . VJZ I VJ /x IVJ 


-2 


HT H9-91 

111 t 1 CJ iJ I 


9.5217 


n K1074-2 0793x10" 


3 


4.91±1.47 


4.36 


-0 00864 

VJ . VJVJOVJrt 


2 26S4+2 2532x10" 

Z . ZVJ*Jte_1_ Z . Z vJtlJZ /x 1 VJ 


-2 


TTpTT ("13 111 

lll_.ll 1 1 ■ ) 111 


9.7068 


q 7011-1-1 9 t ^ c iQvin~ 
y . 1 uii ^ i.zujoy a iu 


■3 


1 fiQ-L-0 ^0 


5.25 


00^1 ^ 
u . uuoio 


1 99Kfi-M 170Qv10~ 
1 . zzou—iz 1 . 1 1 uy A IU 


-2 


TTT ("18 Ql 


9.8470 


y . ooojii . oy i o a i_u 


3 




4.68 


00 

-U . UUdOU 


9 9fi97-l-1 KQ80vlO~ 
z . zvjz 1 _iz 1 . oy OU A IU 


-2 


HI (17-9) 


10.2613 


ll/-^1'J1_Lii.OU10 A 1VJ 


-3 


q q^+1 07 


4.29 


-0 00804 

U . VJVJOUrt 


2 8S80+2 5516x10" 

Z . OvUU 1 Z.OvJIVJ ^x IVJ 


-2 


HI (12-8) /Hell (24-16) 


10.4992 


10. 4871+5. 5642x 10' 


-4 


7.50+0.44 


20.79 


-0.01486 


2.9583+5.5727x10" 


-2 


fSTVl /rOnTTl? 

1 kJl VI / 1 VJUIIJ I 


1 521 


1 5424+4 91 46x 1 


-4 


111 -i-o ^2 

1 • 1 J l_ W . 'Ji 


14.62 


-0 00224 

VJ . VJVJZ Zrr 


1 0418+1 5577x10" 

l.UHO-l L.OtJI / A IU 


-2 




10.6822 


10. 6662+1. 1298x 10 


-3 


6 45+0 83 

VJ • rtrj 1— V ' . < 1 'J 


9.89 


-0 01 386 

VJ . VJ 1 • J'. ' V t 


2.7287+1.1256x10" 


-2 


HT ("16-91 

111 1 1 VJ <J ) 


1 8036 


10 7851+1 8913x10' 

1U. f O'Jl — I L.UJlUA 1 VJ 


-3 


2 "39+0 61 

Z . -J C _l_ V J . VJ 1 


5.44 


-0 00462 

U . VJVJUJZ 


2 2775+1 8887x10" 

Z . Z f 1 ij_l L . OOO f A IU 


-2 


HT ("25- ml 

111 \ ZvJ IVJ / 


10.8514 


10 8588+7 4041 x10' 


-4 


10 21 +0 74 

1VJ.Z1 — l_VJ. 1 rt 


16.24 


-0 01979 


3 1 739+7 4205 X 1 0" 

tJC_L 1 .IZVJvJ/XlVJ 


-2 


HT (23-101 

111 \ ZrJ IVJ J 


11.2429 


11.2206+1.3047x10' 


-3 


1 66+0 24 

1 . V ' \ 1 L_ V ' • Z1 


5.39 


-0 00237 


3 0392+1 4575x10" 

•J . VJvJ C Z _1 L.1U 1 tj A IU 


-2 


HI (9-7)/HeII (23-16) 


11.3026 


11.2895+7.4253x10 


-4 


22.33+1.85 


15.21 


-0.04125 


2.7936+7.4417x10" 


-2 


HI (22-10)/[CaV]? 


11.48200 


11.4801+1.7865X lO- 


-3 


1.59+0.49 


5.40 


-0.00257 


1.8401+1.9194x10" 


-2 


HI (15-9) 


11.5395 


ll. 5197+4.6056X lO- 


-4 


5.37+0.35 


21.94 


-0.00799 


2.2187+4.7369x10" 


-2 


Hell (20-15) 


11.7179 


ll. 7268+1. 5285x10 


-3 


1.85+0.51 


5.98 


-0.00304 


1.7858+1.6108x10" 


-2 


HI (21-10) 


11.7914 


11.7675+1.0107x10 


-3 


3.93+0.59 


10.03 


-0.00605 


2.1741+1.0708x10" 


-2 


HI (20-10) 


12.1568 


12.1338+4.4487x10 


-4 


3.50+0.19 


24.69 


-0.00546 


2.5959+4.6772x10" 


-2 


HI(7-6)/HI(ll-8)/HeII(14-12) 


12.3669 


12.3541+5.9849x10 


-4 


70.95+3.76 


26.25 


-0.12208 


3.5277+6.0072x10" 


-2 


HI (14-9) 


12.5837 


12.5676+1.3871x10 


-3 


14.57+1.72 


9.39 


-0.02618 


3.6829+1.4051x10" 


-2 


HI (19-10) 


12.6110 


12.6248+8.5191x10 


-4 


6.16+0.85 


10.99 


-0.01084 


1.9122+8.4719x10" 


-2 


PAH 


12.7000 


12.6986+1.2823x10" 


-3 


7.76+0.85 


11.18 


-0.01320 


3.8021+1.4296x10" 


-2 


[Nill] 


12.7288 


12.7362+4.4682x10' 


-4 


5.73+0.41 


22.00 


-0.00813 


2.0217+4.7863x10" 


-2 


[Nell] 


12.81355 


12.8281+8.1494x10" 


-4 


2.16+0.32 


10.87 


-0.00398 


1.7358+8.1869x10" 


-2 


HI (18-10) 


13.1880 


13.1901+1.6532x10 


-3 


3.77+1.92 


5.38 


-0.00332 


1.1475+1.9459x10" 


-2 
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TABLE 3 

Fitted parameters of IGR J16318-4848 SOFl-VISIR-Spiteer SED a . 



Model 


T* (K) 


R.» 


- 1 - rim 


hrim 


Tin (K) 


Tin 


^out 


q 


X 2 /dof 






«0 
kpc 


K 


kpc V V wb K 1 )) 


K 


kpc V ( cub \ l l) 


kpc V l LOA W,l 






Flaring disk 


20000 


11.72±0.55 


3791±253 


0.33±0.22 


1716±535 


0.08±0.05 


3.28±0.15 


0.5 


1.20 (1259) 


Viscous disk 


20000 


11.71±0.49 


3786±219 


0.037±0.002 


767±32 


0.74±0.06 


3.47±0.10 


0.75 


1.19 (1259) 


Viscous disk 
Viscous disk 


20000 
20000 


11.20±0.06 
9.86±0.08 


4000 (fixed) 
4500 (fixed) 


0.037±0.002 
0.036±0.001 


790±19 
834±20 


0.70±0.03 
0.64±0.03 


3.49±0.07 
3.50±0.08 


0.75 
0.75 


1.20 (1260) 
1.19 (1260) 


Viscous disk 
Viscous disk 
Viscous disk 


20000 
20000 
20000 


8.16±0.12 
5.80±0.20 
(fixed) 


5000 (fixed) 
5500 (fixed) 
5963±35 


0.035±0.001 
0.035±0.001 
0.035±0.001 


868±21 
895±21 
912±22 


0.60±0.03 
0.57±0.02 
0.56±0.02 


3.51±0.07 
3.53±0.07 
3.55±0.07 


0.75 
0.75 
0.75 


1.20 (1260) 

1.21 (1260) 

1.22 (1260) 


Viscous disk 


20000 





6000 (fixed) 


0.035±0.001 


923±20 


0.55±0.02 


3.55±0.07 


0.75 


1.22 (1260) 


2-sphere 






Tdusti (K) 


Kdust^ 
D„ 


Tdust 2 (K) 


Kdust2 
D» 








component 


20000 


13.25±0.10 


2668±43 


0.32±0.01 


379±3 


3.02±0.02 






1.29 (1260) 



a The best fit with Sofl-VISIR- Spitzer data is obtained for the viscous disk model, for T r i m — 3786 K, as shown in Figurc[T] However, this value 
being low for a sgB[c] star, we also tried various fits with frozen T r i m going from 4000 to 6000 K. The fit, which physically better describes a sgB[c] 
star with T r i m ~ 5500 K, is reported in Figure [8] Since fitting with the 2-sphcrc component model gives a higher \r than fitting with the flaring 
disk or viscous disk models, the Sofl-VISIR- Spitzer data allow us to exclude the spherical repartition of the dust surrounding IGR J16318— 4848. 



TABLE 4 

Fitted parameters of IGR J16318-4848 SOFl-Spitzer SED a . 



Model 


T. (K) 


R» 
D, 


T • 

± rim 




Tin (K) 


'in 


^out 


H 


Xl (dof) 




fixed 


R e 

kpc 


K 




K 


!£vW*)) 




fixed 




Flaring disk 


20000 


12.22±0.38 


3538±180 


0.038±0.006 


629±48 


0.80+/0.17 


3.08±0.35 


0.5 


1.30 (994) 


Viscous disk 


20000 


12.23±0.45 


3532±177 


0.022±0.002 


587±31 


1.37±0.18 


3.29±0.23 


0.75 


1.30 (994) 


Viscous disk 
Viscous disk 


20000 
20000 


11.17±0.06 
9.79±0.08 


4000 (fixed) 
4500 (fixed) 


0.023±0.002 
0.023±0.001 


640±19 
679±18 


1.14±0.08 
1.01±0.06 


3.22±0.20 
3.20±0.18 


0.75 
0.75 


1.31 (995) 

1.32 (995) 


Viscous disk 
Viscous disk 
Viscous disk 


20000 
20000 
20000 


8.05±0.12 
5.59±0.20 
(fixed) 


5000 (fixed) 
5500 (fixed) 
5922±33 


0.023±0.01 
0.023±0.01 
0.023±0.01 


707±19 
729±19 
742±19 


0.93±0.05 
0.88±0.05 
0.85±0.05 


3.19±0.16 
3.19±0.15 
3.19±0.4 


0.75 
0.75 
0.75 


1.33 (995) 

1.34 (995) 

1.35 (995) 


Viscous disk 


20000 





6000 (fixed) 


0.023±0.01 


759±18 


0.82±0.04 


3.19±0.15 


0.75 


1.36 (995) 


2-sphere 






Tdusti (K) 


Kdusti 
D» 


Tdust 2 (K) 


Kdust2 








component 


20000 


13.03±0.17 


3108±84 


0.27±0.01 


412±3 


2.82±0.02 






1.32 (995) 



k The fits with Sofl-Spitzer data have a higher Xr than the fits including Sofl-VISIR and Spitzer data, and do not allow to exclude the 2-sphcrc 
component model. 



